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The explicit construction of low dimensional models with Jastrow ground-state wave- 
functions has attracted considerable recent interest |I|-§]. In one dimension, Shastry and 
Haldane |5||| have demonstrated that the ground-state of the ID Heisenberg model with a 
1/r 2 exchange interaction is a Gutzwiller state for the half filled infinite-?/ Hubbard model. 
Haldane has shown how the spectrum of this model can be written in terms of a generalized 
type of Jastrow wavefunction with excitations of novel statistics . 

Kuramoto and Yokoyama |8| have recently extended these results to include holes, 
demonstrating that the corresponding ID supersymmetric t-J model is also characterized 
by a Gutzwiller ground-state. Most recently, Kawakami has obtained an asymptotic Bethe- 
ansatz (ABA) solution for the model, based on the observation that the ground-state wave- 
function is a product of two-body functions |J. Assuming factorizability, he derived the 
spectrum of the system, which was conjectured to be exact. The low-energy critical be- 
havior of the model has been identified as a Luttinger liquid |TU , II| ; the spin and charge 



excitations are described independently by c = 1 conformal field theories. 

In the case of the 1/r 2 Bose gas JT2], and the Shastry-Haldane 1/r 2 Heisenberg chain [|7|], 



the ABA has been shown to furnish the correct spectrum, despite the long-range nature of 
the interactions. A remarkable feature of these models is that excited states are obtained 
from the ground-state by introducing zeros into the Jastrow wavefunction, in a manner 
reminiscent of Laughlin's description of quasiparticles in the fractional quantum Hall effect. 
This motivates us to examine the 1/r 2 supersymmetric t-J model in a similar vein. Here, 
we show how this philosophy can be used to construct the excited state Jastrow wavefunc- 
tions of the 1/r 2 supersymmetric t-J model and indeed, the spectrum confirms Kawakami's 
conjecture. In addition to the spectrum, we are able to obtain the spin degeneracies of each 
state, permitting us to write the the free energy in closed form. 
The Hamiltonian for the one-dimensional t-J model is given by 

H = [-^laCja] + [ J ij ( S * ' S i - \ n i n o)] . C 1 ) 
where we implicitly project out any double occupancies. We take t^- = = t/d 2 (i—j) where 



d(n) = — sin(wr/iV) is the chord distance consistent with periodic boundary conditions on 
N lattice sites PH . 

States in the Hilbert space can be represented by spin and hole excitations from the 



fully-polarized up-spin state \P) [|14]]. If we let Q denote the number of holes and M denote 
the number of down-spins, then S z is given by S z = (N — Q)/2 — M. The wavefunctions are 
given by 

iv>>=£v^)n^iKj p >> ( 2 ) 

x,y a i 

where the amplitude ip(x, y) is symmetric in x = (xi, X2, ■ ■ ■ , the positions of the down- 
spins, and antisymmetric in y = (2/1 , 2/2 , • • • ,Vq), the positions of the holes. S~ a = cj, iC^t 
is the spin-lowering operator at site x a and h) y . = c Vil creates a hole at site yi. 

We can construct a general class of states corresponding to states of uniform motion 
and spin polarization. To describe these states, we generalize Kuramoto and Yokoyama's 
Jastrow ground-state || as follows 



ipc{x,y, J s , Jh) = exp 



2iti 
1^ 



Js £ X a + Jh £ Vi 



^o(x,y)= Y[ d 2 (x a - x p ) J[d(yi - yj) \{d(x a - y 4 ). (3) 

a</3 i<j a,i 

Here, J s and Jh govern the (uniform) momenta of down-spins and holes respectively. J s and 
Jh take on either integral or half-integral values as appropriate to insure that i\)q has the 
correct periodicities under x a — > x a + iV and yi ^ yi + N. 

The Hamiltonian can be broken up into four parts, H = + H° + H mt , where 

(T*) is the up (down) spin transfer operator, H° is the spin exchange operator and H mt is 
the diagonal interaction term. When H acts on i/)q, only affects the y variables and H° 
only affects the x variables. As a result, these operators are easy to treat and yield only two 
and three body terms when appropriate conditions on J s and Jh are met |5|,|6|||. 

However, because exchanges pairs of x a and yi, this term must be treated differently. 
In general, it is not true that T^ipc) — T^^g) because does not commute with the spin 
raising operator. This difficulty was overlooked in earlier work ||. To deal with T^, we use 



an alternate representation for \i[)q) in terms of up-spins and holes. Let us introduce the 
N — M — Q coordinates u = (ui, u 2 , ■ ■ ■ , Un-m-q) which give the location of the up-spins. 
Wavefunctions in this representation are given by the spin rotated version of Eq. (|2]) where 
the x are replaced by u and M is replaced by N — M — Q. Making this transformation, we 
find 

N 

i/>g(x, V, Js, Jh) = Aip G (u, y;N- J s , J h - J s + y), (4) 

where the set of N coordinates (x, y, u) exhausts the entire lattice. A is a constant indepen- 
dent of the spin and hole coordinates. Using this identity, the down-spin transfer operator 
gives 

ipa(x,y) *pG(u,y) 

and can thus be treated in a similar manner as T*. The result gives two and three body 
terms in the variables u and y. These terms can then be converted into sums over the x and 
y variables by making use of the fact that (x, y, u) runs over the entire lattice. 

When the separate terms that contribute to the Hamiltonian are combined, we find that 
the two body terms drop out and the three body terms combine to give constants. As a 
result, tpG with total momentum P = ^r(J s M + J h Q) is an exact eigenstate of H with energy 



N 

N 2 2 
+Q 



2 E = ^M(M 2 - 1) - 2MJ S (N - J s 



\(N 2 - 1) + \(Q 2 - 1) + \{M + Q)(2M - Q) 



-2J h (N - J h ) + 2(J S - J h 



(6) 



The cancellation of the many body terms, and thus this result, is only valid under the 
conditions \ J S - N/2\ < N/2 — (M — 1 + Q/2), \J h - N/2\ < N/2 - (M + Q - l)/2 and 
\Jh — Js\ < {M + l)/2. For a given S z , the minimum energy is given when J s and Jh 
are as close to N/2 as possible. The ground state is given when S z is either or 1/2 
and is a singlet whenever possible ||. When Q = this reduces to the result for the 
Heisenberg chain |)]||. These energy levels have also been found by Ha and Haldane [15 



where J T = J h - N + (M + Q + l)/2 and J i = J h - J s - (M - l)/2. From these energy 



levels we find the spin and charge velocities to be identical to the previous results [B,15 



To investigate the other excited states of the system, we introduce zeros into the 
wavefunction by premultiplying it with polynomials of X a = exp(2nix a /N) and Yi = 
exp(27u?/j/iV). The wavefunctions thus take the following modified Kalmeyer-Laughlin form 



0: 



4>{x,y) = $ a (X,Y)$ b (Y)* , 



(?) 



where $ s and <f> b are completely symmetric under pairwise interchange of their arguments. 
These states will be termed "fully-polarized spinon states" . Loosely speaking, the polynomi- 
als $ s and <&h can be regarded as spin and charge quasiparticle wavefunctions respectively. 

When the Hamiltonian acts on this wavefunction, once again all three-body terms com- 
bine to give constants. However, in this case, some two-body terms remain and we are left 
with the eigenvalue equation 

N 2 



E$ B $ h = E $ B $ b + H 1 + H 2 + H 3 , 



(8) 



where 



H x = 2$ h 



E^E^^, 



E^E^(^ 

i z i<j II 1 3 



H 2 = 4j2di$ s di$ b 



2^E 

a</3 



X n — Xa 



(9) 



Here we denote W = (X, Y) = (X x . . . X M , Y x . . . Y Q+M ) and = W^d/dW^. In deriving 
this, we have shifted $ s by the ground state configuration, W^ 2 . As a result, E is 
given by Eq. (§) with J s = J h = N/2. We require that (degree $ s | < N/2 - (M - 1 + Q/2), 
(degree <5 s 3> h | < N/2 - (M + Q - l)/2 and (degree $ h \ < (M + l)/2, which is to hold for 
each variable X a or Y^ independently. 



The first term, Hi, does not mix <P S and <&h and has been solved by Sutherland ||12|| . 



However, H 2 mixes <P S and ^ and H 3 does not act symmetrically on <P S . As a result, they 
are harder to deal with. We follow Sutherland and start by choosing the following symmetric 
basis functions: 

WW)=EIIiC 

^;W) = Enf'. (io) 

{Pi} i 

where the quantum numbers {ni, . . . , um+q} and {mi, . . . , mg} are taken to be in increasing 
order and {P M } and {Pi} denote permutations of the indices. These quantum numbers are 
integral or half-integral as required by periodic boundary conditions. 

In this basis, labeled by the two sets of quantum numbers {n^} and {nij}, the Hamilto- 
nian, considered as a matrix, can be shown to be upper-triangular. Eigenvalues are found 
by reading the diagonal-elements labeled in terms of the quantum numbers {n^} and {mj 
||17|| . The result simplifies when written in terms of a conjugate set of quantum numbers 
{ Ji, J 2 , . . . , J m +q} and {h, I 2 , . . . , Iq} defined by 



J M = n M + nJ raJ = -(2^-(M + Q)-l) 



/. = m . + m o m o = l(2i_Q_l) > (11) 

where {n^} and {mj must satisfy the conditions specified before. This translates into the 
conditions |J„| < (JV - M + l)/2 and < (M + Q)/2. The energy is 

where the pseudomomenta, p^, are given by the following equations: 

Q m+q 
p^N = 2?r J M - ttJ2 sgn(> M - ft) + vr J] sgn(p M - p„) 
i=i i/=i 

M+Q 

27r/i = 7T J] sgn(ft - p M ). (13) 

The above equations correspond to the asymptotic Bethe-ansatz equations obtained by 
Kawakami ||. Our result thus confirms that the ABA spectrum is exact. 

6 



Here the resulting and {q{\ must lie between — ir and 7r. The set of M + Q distinct 
quantum numbers is in ascending order and governs the spin excitations. We restrict 
them to take values in the range [—(N — M — l)/2, (N — M — l)/2] to guarantee that they 
generate fully spin-polarized states. There are N — M values in this range, of which M + Q 
are occupied and 2S Z are empty. A spin configuration can be represented by a sequence of 
N — M digits such as {S} = (0111001011) s , where 1 represents an occupied quantum number 
and an unoccupied quantum number. These empty values are identified as spinons 0; 
a sequence of 2j r consecutive zeros corresponds to a symmetric bound-state of 2j r spinons, 
thereby creating an excitation of spin j r with spin degeneracy 2j r + 1. On these physical 
grounds, we anticipate a spin degeneracy in the thermodynamic limit given by 

w s = II(2J + l) n0) (14) 

3 

where n(j) is the number of sequences of zeros of length 2j. The set of Q distinct quantum 
numbers I iy in ascending order and taking values in the range [— (M + Q)/2, (M + Q)/2], 
governs charge excitations. 

To complete the study of the model and confirm our interpretation of the quasiparticle 
degeneracies, we looked at exact diagonalization of small systems (N < 10 with holes). As 
an example, the low-lying states of the N = 10, Q = 2 model is shown in Fig. p]. We 
summarize the numerical result as follows: 

1. The spectrum described in terms of the real pseudomomenta {p^} and {qi\ span the 
full set of energy levels of the system. 

2. The real pseudo momentum states are all highest weight states when {p^} 7^ ±7r. 

3. The spin degeneracy rule is obeyed for all internal sequences of zeros. 

Certain small corrections to the spin degeneracy rule apply when there are zeros at either 
end of {S} which we shall not enumerate here, and which are not important in the thermo- 
dynamic limit P~7[] . 
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Finally, we may use the spectrum generated by the "fully-polarized spinon states" and 
the supermultiplicity rule to obtain the free energy of the model in the thermodynamic 
limit. Besides the "particle-state" solutions of equations fll2"|) and (|I5D , we have to take into 
account the "hole-state" solutions ]n|. At thermal equilibrium, the distribution functions 
of these solutions are determined by minimizing the free energy functional fK||2"0| ], F = 
E — TS — fi(N — Q), with the constraint that each "fully-polarized spinon state" described 
by quantum numbers {J^\ and {ij} is associated with a spin degeneracy ws as given in 



([14]), where fi is the chemical potential. Following the standard methods of Takahashi [19 



minimizing the free energy for a given quantum number distribution gives the following free 
energy 



F{T)/N 



T 

-a / dp In 



1 + e 



-/Mp) 



(15) 



where e s is determined by the coupled equations 



2t s (p) = e (p) - 2a - Tin 



1 + e 



-/Mp) 



e c (q) = 2a — Tin 



1 + e 



-/M?) 



(16) 



Here eo(p) = \t{p 2 — ^) + /i and a = |t7r 2 + In the limit of half filling, ji — > oo, 
e s — > t(p 2 — 7r 2 )/4, and the free energy reverts to the form obtained by Haldane for the 
corresponding Heisenberg model [fTj^Tj. For general /z, elimination of e c (p) yields the result 



e s (p) = e (p) - Tin 



I + 2e^( e o(p)+ a ) cosh(/5a; 



1/2' 



(17) 



We have verified that high temperature expansion of this free energy in powers of (3 correctly 
reproduces the first two non-trivial terms in the high temperature perturbation theory. 

From the free energy, it is not clear whether we may make a unique identification of the 
statistics of the spin and charge excitations. We note that the 5 = 1/2 spinon excitations 
always combine into a state with a symmetric spin wavefunction, thus 2S spinons form a 
state with total spin S. In the limit of zero doping, the free energy is that of spinless fermions 
|7|j2l[|. We can equally well regard the spinon excitations as S = 1/2 fermions in a state with 



S 



a fully antisymmetric spatial wavefunction; or alternatively, as hardcore S — 1/2 bosons in 
a fully symmetric spatial state. 

In summary, we have derived the spectrum of the ID t-J model with 1/r 2 long-range 
exchange and hopping by the introduction of zeros into Jastrow ground-state wavefunctions. 
Our solution confirms Kawakami's conjecture that the ABA provides the exact spectrum, 
suggesting that despite the long-range nature of the interactions, two-body scattering domi- 
nates the long-wavelength physics. By interpreting multiple occupancy of momentum states 
in the spinon wavefunction as symmetric bound-complexes of spinons, we have been able to 
determine the degeneracies of the states needed to construct the free energy. Further work 
is required to determine the integrability conditions of this model. There are also several 
possible generalizations: most notably, SU(N) generalizations and the appealing possibility 
of Jastrow-integrable impurity models. 
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work was supported in part by the U.S. Department of Energy under Grant No. DE-FG05- 
85ER-40219, the National Science Foundation under Grant NSF-DMR-13692, and by a Sloan 
Foundation Grant. 
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FIGURES 

FIG. 1. Low- lying energy levels of the 10 site 2 hole system from exact diagonalization. The 
numbers associated with each state list the spin degeneracies starting with spin on the left. For 
example, the number "331" indicates that we have 3 states with S = 0, 3 states with S = 1 and 1 
state with 5 = 2. 
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